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FOREWORD 


This  report  was  prepared  under  Contract  AF  0H(635)-Z  1  55,  "Study 
of  Target  Penetration  Prediction  By  High  Speed  and  Ultra  Higli  Speed 
Ballistic  Impact",  AFSC  Project  ‘hS60.  Work  uas  administered  under 
the  direction  of  APGC  (PGWR  i  ),  EglinAFB,  Florida,  with  Mr.  A.  G. 
BiU'k  as  Proj(“ct  Engineer. 


Catalog  t\'irf/,S'  with  nbstrncts:  rn/iy  hi'  ioiind  in  thi'  h;ii-k  o!  this  documrnt. 


Complete  results  of  the  preliminary  statistical  analysis,  which 
were  partially  reported  in  the  previous  quarterly  report,  are  reported 
herein.  This  statistical  correlation,  based  on  all  hypervelocity 
terminal  ballistic  data  gathered  prior  to  December  1961,  attempts  to 
relate  the  depth  of  penetration  in  semi-infinite  targets  with  ten  inde¬ 
pendent  variables.  Without  any  initial  assumptions  being  made  regarding 
the  process  of  ballistic  impact  or  the  shape  of  the  craters  formed,  this 
analysis  produced  an  equation  quite  similar  to  the  empirical  equation 
used  by  many  Investigators  to  fit  their  data. 

Also  reported  are  initial  attempts  to  formulate  a  theoretical  model 
for  the  purpose  of  testing  accumulated  experimental  data.  This  theoret¬ 
ical  model  is  developed  from  a  consideration  of  energy  conversion  during 
impact. 
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LIST  OF  SYMBaS 

P  crater  depth  measured  from  original  target  surface 
c 

D.  crater  diameter  at  original  target  surface 
Vg  crater  volume  below  original  target  surface 

Dp  projectile  diameter 

p  density 

V  impact  velocity 

c  dilatational  wave  velocity  in  a  given  material 

m  mass 

T  target  temperature 

V  volume 

p  pressure 

A  area 

t  time 

E  energy 

€  strain 

V  yield  point 

U  ultimate  strength  (tensile) 

L  latent  heat  of  fusion 

SUBSCRIPTS 

t  target 

p  projectile 

c  crater 

o  SIP  conditions  or  arbitrary  reference  point 
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The  purpose  of  this  study  is  to  gather  and  assemble  existing  data 
on  ballistic  impact  and  on  material  failure,  especially  at  high  impact 
velocities  or  large  loading  -  to  establish  the  relative  importance  of 
such  factors  as  projectile  velocity,  mass,  sectional  density,  projec¬ 
tile/target  contact  areas,  etc.,  and  to  use  this  information  to  deduce 
the  mathematical  relationships  of  critical  factors  as  the  target  struc¬ 
ture  responds  to  impact  and  is  penetrated. 

Existing  experimental  data  relative  to  ballistic  Impact  at  high 
velocities  are  being  evaluated  on  a  statistical  basis  through  the  use 
of  an  RPC  4000  digital  computer.  The  general  form  of  the  statistical 
approach  was  outlined  in  the  First  Quarterly  Report^.  A  preliminary 
analysis  of  the  correlation  between  depth  of  penetration  and  ten  inde¬ 
pendent  variables  was  reported  in  the  Second  Quarterly  Report^  and  will 
be  discussed  more  fully  in  this  report. 

In  addition,  the  general  areas  of  ballistic  impact  and  material 
failure  are  being  investigated  in  order  to  develop  relationships  which 
may  be  tested  against  existing  experimental  data.  Some  general  aspects 
of  target  behavior  under  ballistic  impact  were  discussed  in  the  Second 
Quarterly  Report  and  some  justification  was  given  for  the  use  of  static 
or  quasi-static  material  parameters  in  the  initial  statistical  analysis. 


EMPIRICAL  MODEL 

The  initial  approach  to  the  empirical  model  was  a  statistical 
analysis  of  the  factors  influencing  the  depth  of  penetration  (P(.)  in 
semi-infinite*  targets.  Penetration  depth  was  considered  a  function 
of  ten  independent  variables  according  to  the  equation 

Pc  =  ko  v  *^1  /Ot  *'2  V  *^4  T  *^5  Y.  *^6  Y  ^7  C  u  *^9  D  *^10 

CO  /t/-^p  p  t  p^^p 

This  equation  was  reduced  to  its  associated  linear  form 

*  Semi-infinite  according  to  the  rule-of-thumb  published  by  Kinard,  et  al^, 
that  the  penetration  is  no  greater  than  20X  of  the  target  thickness. 
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log  Pg  =  log  kg  +  log  v  + 


+  k^o  log  Dp 


and  the  method  of  least  squares  was  used  to  determine  a  "best"  set  of 
values  for  the  variable  exponents  (coefficients)  kj  based  on  minimizing 
the  value  of 


lO 


^  !  I09  ^^c ^experimental  ^^c^calculated 
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1--0  J 

The  independent  variables  were  dropped  one-by-one  starting  with  Dp 
the  effect  on  the  kj  is  summarized  in  Table  I* 


and 


After  the  k^  were  computed,  various  statistical  tests  were  made  to 
determine  the  validity  of  the  results  obtained.  The  per  cent  of  the 
variance  in  the  log  of  dependent  variable  (P^)  that  is  explained  by  the 
logs  of  the  independent  variables  is  given  by  the  multiple  correlation 
coefficients  listed  in  Table  I.  The  significance  of  these  multiple 
correlation  coefficients  was  discussed  in  the  previous  progress  report^. 
An  analysis  of  this  variance  using  the  F  distribution  tests  the  null 
hypothesis  that  all  /3  j  are  zero.  Separate  tests  of  the  individual 
hypothesis  /3<  =  0  were  made  using  the  t  distribution.  Variables 
corresponding  to  the  /3^  that  were  not  significant  by  the  preceding 
test  were  then  removed  and  separate  tests  were  made  on  each  of  the 
removed  variables  to  see  whether  the  adding  of  this  variable  improved 
the  fit.  Individual  correlation  coefficients  between  all  independent 
variables  were  also  calculated. 


The  significance  of  the  target  yield  strength  above  many  of  the 
other  variables  contributing  to  the  multiple  correlation  coefficients 
was  noted  in  the  previous  progress  report.  However,  the  target  yield 
strength  (Y^),  the  dilatational  wave  velocity  in  the  target  (c^),  and 
the  target  ultimate  strength  (U^)  are  statistically  indistinguishable  as 
evidenced  by  the  fact  that  the  calculated  individual  correlation  between 
them  is  over  95%.  The  diameter  of  the  projectile  (D  )  is  also  statis¬ 
tically  indistinguishable  from  the  volume  of  the  projectile  (Vp)  so  that 
the  effect  of  dropping  (Dp)  should  be  insignificant.  With  the^exception 
of  the  above  correlations  and  the  dependence  of  on  the  rest  of 
the  individual  correlation  coefficients  show  a  low  degree  of  interdependence. 


The  close  correlation  between  Dp  and  V  and  between  U^,  and 
confuses  the  kj  calculated  in  the  first  three  lines  of  Table  I.  It  is 
interesting  to  note  that  that  the  last  line  of  the  table  yields 


P 


c 
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with  being  statistically  insignificant  according  to  the  t  distribution. 
This  is  in  good  agreement  with  the  many  investigators'^’^’®*^’®’^’^^’^^ 
who  found  empirically,  that  the  volume  of  the  crater  is  proportional  to 
the  projectile  kinetic  energy 


Vc  --  k  />p  \  v") 


(2) 


coupled  with  the  observation^*®’^^’^^’^^’^^  that  the  craters  formed  are 
hemispherical 


V 


c 


(3) 


since  these  two  observations  yield  an  equation  of  the  general  form 


2  1  1 

P  )c  v'  3  /D  '  "3  u  3 
C  /  p  '^p 


(4) 


The  agreement  in  the  exponents  of  v  and  V  between  equations  (l)  and  (4) 
is  quite  good.  However,  the  disagreement^between  the  exponent  of  is 
quite  substantial.  Palmer^^  points  out  that  all  Utah  cratering  '  P 
studies  seem  to  best  fit  an  equation  similar  to  equation  (4)  in  which  the 
exponent  of  is  but  this  does  little  toward  resolving  the  disagree' 
ment.  Since  the  multiple  correlation  coefficient  shows  that  87^  of  the 
variance  in  (P(.)  may  be  explained  by  the  variance  in  v,  /2p,  and  Vp,  it 
is  not  surprising  that  equation  {4)  has  been  quite  successful  in  fitting 
individual  experimental  data. 
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Palmer'*’  observes  that  the  empirical  crater  depth  seems  to  be  in¬ 
versely  proportional  to  the  _i_  power  of  the  target  shear  strength  and 
Eichelberger^^  observes  that  ^it  seems  to  be  inversely  proportional  to 
the  power  of  the  target  Brinell  hardness.  Both  the  Brinell  hardness 
and  the  shear  strength  are  quite  closely  correlated  with  the  yield  strength, 
and  5  of  Table  I  indicate  a  dependence  of  Pq  on  either  the  plus 


and  lines  4 

or  minus  -X 
3 


power  of  the  target  yield  strength. 


The  preliminary  analysis  was  made  on  approximately  1200  experimental 
shots  only  half  of  which  were  reported  in  enough  detail  to  enable  assign¬ 
ment  of  strength  parameters  to  the  target  and  projectile.  Of  the  resulting 
600  shots,  Only  174  were  performed  at  velocities  greater  than  the  bulk  wave 
velocity  iri  the  target  medium  and  over  90^  of  these  "supersonic"  shots  were 
in  lead  targets  where  the  bulk  wave  velocity  is  quite  low  (about  5400  feet 
per  second).  This  analysis,  then,  seems  to  indicate  that  the  crater  depth 
can  b'=»  adequately  predicted  on  the  basis  of  4  or  5  out  of  the  ten  indepen¬ 
dent  variables  chosen,  at  least  for  impact  velocities  below  the  bulk  or 
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dilational  wave  velocity  in  the  target  medium. 

A  second  statistical  analysis  was  attempted  during  this  report 
period  using  only  the  174  shots  at  impact  velocities  above  the  bulk 
wave  velocity  in  the  target  medium.  The  results  of  this  analysis  were 
completely  inconclusive  due  primarily  to  the  fact  that  all  but  17  of 
the  shots  were  fired  into  lead  targets.  This  overwhelming  number  of 
lead  shots  created  additional  correlation  between  the  various  material 
parameters  yielding  a  spurious  set  of  results.  Since  this  analysis  was 
performed,  the  results  of  about  60  shots  using  aluminum  projectiles  on 
aluminum  targets  at  impact  velocities  between  27,000  and  33,000  feet/ 
second  have  been  received.  A  new  high  impact  velocity  analysis  will  be 
initiated  soon,  subject  to  the  boundary  conditions  outlined  below. 

Our  empirical  analysis  to  date  clearly  indicates  that  we  were  too 
concerned  with  the  quantity  of  experimental  data  and  not  enough  concerned 
with  quality.  An  analysis  of  (Pc/Op)  as  outlined  in  the  First  Quarterly 
Report  is  almost  completed  on  all  snots  below  the  bulk  wave  velocity  in 
the  target  medium.  Before  any  further  analyses  are  started,  the  terminal 
ballistic  data  will  be  very  carefully  screened  according  to  the  conditions 
under  which  it  was  taken  and  care  will  be  exercised  in  balancing  the  data 
to  some  degree  so  that  a  more  even  distribution  of  target  and  projectile 
materials  and  Impact  velocities  will  be  used. 

Spurious  results  caused  by  tne  statistical  interdependence  of  inde¬ 
pendent  variables  clearly  illustrates  the  need  of  a  theoretical  analysis 
to  complement  any  empirical  study.  The  theoretical  model  which  follows 
is  merely  an  attempt  to  develop  a  logical  set  of  independent  variables 
in  dimensionally  correct  groupings  against  which  to  test  the  experimental 
data. 


THEORETICAL  MODEL 

In  conjunction  with  a  broad  statistical  study  of  medium  arJ  high  speed 
ballistic  impact  data,  theoretical  studies  of  target  response  are  also  in 
progress.  The  study  is  attempting  to  correlate  target  behavior  with  per¬ 
tinent  impact  and  material  parameters  through  fundamental  considerations 
of  both  the  energy  and  momentum  processes  involved. 

The  mechanics  of  impact  behavior  involves  a  combination  of  energy 
absorbing  and  momentum  absorbing  processes.  In  the  case  of  lower  velocity 
impact,  the  target  tends  to  steadily  resist  the  force  of  impact,  resultin'^ 
primarily  in  an  energy  absorbing  process,  as  the  empirical  models  seem  to 
indicate  (V^.  proportional  to  v^).  This  is  an  over-simplification,  of  course, 
but  the  same  behavior  would  be  expected  to  occur  at  higher  impact  velocities, 
although  to  a  lesser  degree. 


In  higher  speed  impacts  the  ejection  of  target  particles  is  more 
and  more  in  evidence,  and  it  seems  apparent  that  the  high  momentum  of 
the  backsplash  material  indicates  the  importance  of  considering  the 
momentum  absorbing  processes  involved.  This  would  seem  to  imply  that 
the  crater  formation  in  ultra-high  speed  impact  tends  to  become  a  closer 
function  of  velocity  to  the  first  power  as  the  theories  of  Bjork  and  (^ik 
indicate. 

The  kinetic  energy  of  a  projectile  impacting  a  semi-infinite  target 
is  transformed  through  a  number  of  complex  processes.  For  the  purposes 
of  obtaining  approximate  expressions  for  the  principle  energy  absorbing 
processes,  the  projectile  kinetic  energy  can  be  assumed  to  be  absorbed 
as  followsi 

1.  Kinetic  energy  of  the  material  backsplashed  during  cavity 
formation 


Eg  -  (i  mv2)  backsplash 


Various  opinions  have  been  expressed  as  to  the  relative  magnitude 
of  the  kinetic  energy  of  this  ejected  material.  It  does  appear 
that  the  mass  of  the  backsplash  material  is  usually  much  greater 
than  that  of  the  original  projectile.  Although  the  Utah  group^^ 
has  measured  very  high  velocities  (up  to  20  times  the  impact 
velocity  for  micron-sized  spray  particles),  the  average  velocity 
of  the  backsplashed  material  must  be  quite  small.  This  is 
verified  by  photographs  of  impacts  which  show  large  chunks 
of  material  moving  backward  from  the  point  of  impact  with  very 
small  velocities. 

2.  Energy  required  to  melt,  vaporize  or  otherwise  dissociate  projec¬ 
tile  and  target  material  during  the  crater  formation.  It  is 
considered  of  little  consequence  here  that  the  actual  :itechanism 
of  material  separation  is  at  least  in  part  a  result  of  cavita¬ 
tion  following  the  disappearance  of  the  projectile.  Bromberg 
has  pointed  out  that  the  terminal  ballistic  data  of  Summers  and 
Charters^  can  be  explained  using  the  heat  of  fusion. 


=  V, /Ot  Lt  +  Vp  ^pLp  (6) 


If  any  vaporization  takes  place,  it  could  account  for  a 
very  large  portion  of  the  impact  kinetic  energy.  A  simple  cal¬ 
culation  for  an  aluminum  projectile  striking  an  aluminum  target 
at  an  impact  velocity  of  9  km. /sec.  shows  that  it  would  take  l/S 
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of  the  impact  kinetic  energy  to  vaporize  the  projectile. 

If  the  crater  volume  were  57  times  the  projectile  volume 
(not  an  unreasonable  figure)  the  energy  necessary  to 
vaporize  the  entire  crater  would  be  over  10  times  the 
impact  kinetic  energy.  The  vaporization  term  has  been 
eliminated  from  equation  (6)  because  it  appears  that  only 
a  small  portion  of  the  projectile  is  actually  vaporized. 

To  create  the  spray  of  micron-size,  backsplashed  particles 
being  studied  by  the  Utah^^  group  would  require  at  least  a 
portion  of  the  vaporization  energy,  but  here  again,  the  total 
mass  of  this  material  is  small. 

3.  Energy  emitted  as  electromagnetic  radiation  at  all  wavelengths 
due  to  the  heat  generated  at  the  point  of  impact.  This  would 
be  given  by  the  Stefan-Boltzmann  Law  as 


Er  =  e  (T  T'^  At 


(7) 


where  (e)  is  the  emissivity  of  the  target  material,  (  CT  )  is 
the  Stefan-Boltzmann  constant,  (T)  is  the  Kelvin  temperature, 
(a)  is  the  area  of  the  emitting  surface,  and  (t)  is  the  time 
over  which  this  emission  takes  place.  Picking  a  good  value  for 
(T)  is  quite  difficult  although,  if  (a  )  is  assumed  to  be  the 
crater  area,  (T)  might  be  approximately  the  melting  temperature 
of  the  target  material. 

Since  (T)  and  (a)  will  both  be  functions  of  time,  it  might 
be  well  to  define  as 


Er  =  e 

where  (tf)  is  the  total  time  of  crater  development,  and  try  to 
find  reasonable  values  for  the  variations  of  (T)  and  (a)  with 
time.  In  any  case,  this  term  will  probably  be  small  in  compari¬ 
son  with  some  of  the  other  terms. 

4.  Energy  transferred  to  atomic,  molecular,  or  granular  change  such 
as  the  energy  of  recrystallization 


=  Uiuirtc  riniiU'd  quantity 


T(t)''  A(t}  dt  (8) 


15 

The  Utah  group  has  shown  that  from  12  to  15^  of  the  energy 
required  to  push  a  lead  ball  into  a  lead  target  quasi-statically 
goes  into  recrystallization.  It  has  not  yet  been  determined  what 
independent  variables  might  enter  into  this  term. 
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strain  energy  (E^)  absorbed  by  the  solid  target  material. 

This  can  be  approximated  by  an  integration  involving  the 
product  of  the  strain  energy  absorbed  per  unit  volume  of 
target  material  and  the  volume  of  target  material  so  stressed. 


p  6"  dV 


An  estimate  will  be  required  of  the  portion  of  the  target 
that  is  stressed  and  to  what  extent  it  is  stressed  in  finally 
absorbing  a  portion  of  the  terminal  impact  energy.  This  will 
require  an  understanding  of  the  magnitudes  of  the  pressures 
developed,  the  pressure-time  variations  in  the  target,  and 
the  pressure-deformation  characteristics  (equations  of  state) 
of  the  target  material. 

The  volume  integral  in  equation  (lO)  may  be  transformed 
to  a  time  Integral  by  reference  to  Figure  1. 


dV  4  TT  (ct)  d  (ct) 


Since  the  pressure  (p)  is  a  function  of  time  and  the  strain  (£) 
is  a  function  of  pressure,  equation  (lO)  becomes 


Eg  =  27t  p  (t)  e  (p)  (ct)"  d  (ct)  (12) 


The  dilatational  wave  velocity  in  the  unstressed  target 
material  (c)  is  a  constant.  Although  later  stages  of  the 
pressure  wave  may  travel  through  the  target  material  at  a 
greater  velocity  than  (c)  due  to  the  previous  compression 
of  the  target  by  the  wave  front,  the  wave  front  itself  cannot 
move  faster  than  (c)  through  the  previously  unstressed  target 
material.  This  situation  seems  somewhat  analagous  to  "choked 
flow"  of  fluids.  Equation  (l2)  may  then  be  written 


Es  =  27y  p  (t)  e (p)  0^  t^ 


The  total  time  (t)  during  which  stress  energy  acts  as  a 
means  of  removing  energy  from  the  point  of  impact  is  somewhat 
in^ doubt.  However,  it  seems  most  reasonable  to  assume  that 
(t)  is  the  time  after  impact  at  which  the  projectile  ceases 
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FIGURE  I 
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be  a  causitive  force  in  the  growing  crater.  According  to 
Gehring^^  this  time  will  be  given  by 


Before  the  integral  of  equation  (13)  can  be  evaluated, 
it  remains  to  describe  the  variation  of  pressure  with  time 
and  the  variation  of  strain  with  pressure  over  the  interval 
t  .  Attempts  to  arrive  at  a  reasonable  expression  for  the 
pressure  variation  are  underway.  Once  the  pressure  variation 
is  known,  the  variation  in  the  strain  may  be  evaluated  from 
Hugoniots  similar  to  those  given  by  Bjork^^.  It  appears  that 
regardless  of  the  pressure  variation,  the  integrated  form  of 
equation  (13)  will  be 


(c  t')^  pj^  (15) 

K 

where  (pq)  is  the  maximum  pressure  applied  to  the  target 
material,  (^q)  is  the  ultimate  strain  of  the  material  at 
pressure  (pq),  and  K  is  a  constant  dependent  upon  the  mode  of 
variation  of  p(t)  and  £(p). 

It  should  be  pointed  out  here,  that  a  portion  of  the 
energy  which  goes  into  elastic-plastic  strain  energy  may 
account  for  part  of  the  backsplashed  kinetic  energy  (equa¬ 
tion  5).  This  would  be  due  to  the  "decompression  wave" 
tending  to  push  the  bottom  of  the  crater  backward  after 
the  forward  pressure  has  dropped  essentially  to  zero. 

Assuming  that  all  modes  of  energy  dissipation  have  been  accounted  for 
in  the  above  five  expressions,  conservation  of  energy  would  require  that 


E  + 
L 


E  +  E_ 
A  S 


(16) 


If  the  crater  volume  is  assumed  to  be  proportional  to  the  kinetic  energy 
of  the  projectile,  and  the  shape  of  the  crater  is  assumed  to  be  hemispher¬ 
ical,  an  expression  for  P^,  will  result.  As  soon  as  more  thought  is  given 
to  some  of  these  energy  terms  and  as  soon  as  the  experimental  data  is 
arranged  in  a  better  form  statistically,  the  theoretical  expression  will 
be  fitted  to  the  empirical  data  in  order  to  evaluate  the  constants  in  the 
expression  and  to  find  the  degree  to  which  the  variance  in  or  may  be 
explained. 


10 


The  principle  of  impulse  and  momentum  requires  that  the  impulse 
of  a  force  system  on  a  system  of  particles  during  a  time  interval  be 
equal  to  the  change  in  momentum  of  the  system  of  particles  during  the 
time  interval.  If  the  time  interval  in  question  is  taken  to  be  the 
time  during  which  the  projectile  remains  as  a  causitive  force  (t^), 
then  the  impulse  of  the  resulting  impulse  developed  in  the  target 
plus  any  external  reaction  (R)  developed  during  the  time  interval  can 
be  equated  to  the  difference  in  the  momentum  of  the  ejected  particles 
and  the  projectile  impact  momentum. 


F  At  =  (m  V  -  (m 

t'  .t' 

p  A  dt  -Jr  dt  =  (-  ra  v)backsplash  "  Vp 


r\ 

or  since  (Figure  l)  A  =  ^  4  TT  (ct)'^ 

[qi^p  (ct)"  dt  +  Tr  dt  =  (m  ^)i,,ckspUsh  +  '"p  •'^p 


Assuming  there  is  no  external  reaction  at  the  target  (R  =  O),  a 
ballistic  pendulum-type  experiment  to  measure  backsplashed  momentum  might 
enable  the  variation  of  pressure  with  time  to  be  evaluated  for  insertion 
into  the  energy  equation  (13). 
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